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On Orbital Period Changes in Nova Outbursts 
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ABSTRACT 

We propose a new mechanism that produces an orbital period change during a 
nova outburst. When the ejected material carries away the specific angular momentum 
of the white dwarf, the orbital period increases. A magnetic field on the surface of the 
secondary star forces a fraction of the ejected material to corotate with the star, and 
hence the binary system. The ejected material thus takes angular momentum from 
the binary orbit and the orbital period decreases. We show that for sufficiently strong 
magnetic fields on the surface of the secondary star, the total change to the orbital 
period could even be negative during a nova outburst, contrary to previous expecta- 
tions. Accurate determinations of pre- and post-outburst orbital periods of recurrent 
nova systems could test the new mechanism, in addition to providing meaningful con- 
straints on otherwise difficult to measure physical quantities. We apply our mechanism 
to outbursts of the recurrent nova U Sco. 
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1 INTRODUCTION 

Classical novae are binary systems in which mass is trans- 
ferred from a main-sequence star on to a white dwarf by 
Roche-lobe overflow. The critical amount of mass that can 
be accreted on to the surface of the white dwarf prior 
to an outburst i s a strongly decreas ing function of the 
white dwarf mass man fc Lividll986ft . At this mass limit, 
the temperature and density at the base of the accreted 
layer are high enough for hydrogen to ignite. The tem- 
perature then rises rapidly in a thermonuclear runaway 
jStarrfield. Sparks fc Shavh1 fl988) and the pressure at the 
base of the accreted layer becomes high enough that the ac- 
creted mass (and somet imes a little more) is ejected (e.g. 
iKovetz fc Prialnild 11985ft . 

R ecurrent novae show outbursts at intervals of 10 — 
80 yr (|Warnedfl995l : IWebbink et all 11987ft . To account for 
the short timescale between the outbursts, the white dwarf 
in a recurrent nova system m ust have a mass, Mi, close t o 
the Chandrasekhar limit (e.g. iKato fc Hachlsulll988l . Il989ft . 
We consider in more detail the recurrent nova U Sco that 
most recently erupted in 2010. The evolved companion in 
U Sco has a mass of Mi — 0.88 Mq so the mass ratio, q = 
Mi /Mi = 0.64, is relatively large. 

In Section [2] we begin by re-examining all the previ- 
ously proposed sources of orbital period change during a 
nova outburst. We first describe a simple model where the 



material carries away its specific angular momentum, then 
we include mass accretion on to the companion and frictional 
angular momentum losses as the binary moves through the 
common envelope. In Section [3] we propose a new mecha- 
nism for orbital period change involving the magnetic field 
on the secondary star. In Section [4] we apply our model to 
the outbursts in the recurrent nova U Sco. 



2 OUTBURST MODEL 

We consider first a simple model of the outburst where the 
ejected material carries away the specific angular momen- 
tum of the white dwarf. The non-degenerate mass accu- 
mulated on to the surface of the white dwarf is very thin. 
The envelope is ejected when the pressure at the surface 
of the white dwarf reaches a critical value o f the order of 
Pcrit = 10 20 dyncm- 2 (e.g. lFuiimotolll982al lbl: iMacDonaldl 
1983). The critical amount of mass that accumulates before 
a nova outburst is of order 
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(1) 



GMi ' 

for a white dwarf of mass Mi and radius Ri . Both theory and 
observations suggest that all the the material that has been 
accreted since the last outburst is ejected in the outburst 

The angular momentum of the binary star system is 
given by 
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where a is the separation of the two stars and the angular 
velocity, fi, is given by Kepler's law 

n4 



(3) 



Here the mass of the companion star is M2 and the total 
mass of the system is M = Mi + M 2 . 

If the ejected mass carries away its specific angular mo- 
mentum and all the mass is lost from the system, then the 
angular momentum removed from the system is 



AJ S pec = —Ami aiQ, 



(4) 



where ai is the distance of Mi to the centre of mass of the 
binary 

M 2 



ai 



M 



(5) 



We neglect the spin of the white dwarf because it is not 
coupled to the binary orbit. By differentiating equations ([2]) 
and (|3]) we find the change in angular momentum of the 
binary 
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where Aa is the corresponding change in the separation of 
the system (due to mass lost) during the outburst, AMi is 
the change in mass of the white dwarf, AM2 is the change 
in the mass of the secondary and AM is the change to the 
total mass of the system. For the case where all the ejected 
material is lost from the system, we choose AMi = —Ami 
and AM2 = 0. We take the angular momentum loss from 
the system to be equal to that carried away by the ejected 



mass so that A J : 
we find 



AJspcc. With equations ©, © and (© 
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As mass is lost from the system in the outburst the separa- 
tion increases. 

By differentiating equation Q we find the period 
change during the outburst to be 

ap _ aq 
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and, for the case where the material carries away all of its 
specific angular momentum, with equation JTJl we find 



AP _ Ami 
P M ' 



(9) 



Since Ami > we see that the period of the system should 
increase during the outburst if the material carries away its 
specific angular momentum. In Fig. [T] we plot ^^J M as 
a function of the mass ratio, q, for the different outburst 
models we consider. The solid constant line corresponds to 
this model where the mass ejected carries away its specific 
angular momentum. 



2.1 Mass Accretion on to the Companion 

Now we consider the effect of a fraction of the ejected mass, 
f3, that may be captured by the companion in the outburst. 
We take A Mi = -Ami and AM 2 = 0Ami so that 
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Figure 1. The change 



as a function of mass ratio, 



q = M2/M1. For the solid line the ejected mass just carries 
away its specific angular momentum (equation [S). The dotted 
line shows the case where accretion onto the secondary is consid- 
ered (corresponding separation change in equation [TTJ. The short 
dashed line includes accretion on to the secondary and frictional 
angular momentum losses (corresponding separation change in 
equation 1141 . The long dashed line takes into account the angu- 
lar momentum change of the ejected material due to a magnetic 
field on the secondary that has _Ra = 0.75 a (corresponding sep- 
aration change in equation 1291 1 . We calculate the period change 
from the separation change with equation JSJl. 



The corresponding change to the separation is then 



Aa Ami 
~ = M 



1 + P [2q 
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jShara et all Il986h . where q = M 2 /Ml. This reduces to 
equation ([7]) with j3 — 0. In the absence of strong magnetic 
effects the maximum value of P is the fractional area of the 
companion's radius. We can estimate the captured fraction 
of mass with 
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Because the secondary fills its Roche lobe, we can estimate 
the stellar radius with R2 = -Rl, where 



Rl 



0.49 q- 



(13) 



v 0.6q3 +hi(l+a3) / 

|Eggleton|[l983l ). We compute the period change for this 
model including mass accretion on to the secondary star 
with equations and (jlll) . In Fig. [T] we plot ^_IJ M as 
a function of the mass ratio for this model (dotted line). 
Generally, the separation change is not so large as the pre- 
vious case where all the ejected mass left the binary system. 
However, for the largest mass ratios the period increases 
compared to the previous model. 
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2.2 Frictional Angular Momentum Losses 

Angular momentum is lost from the binary system be- 
cause of frictional angular momentum losses as the binary 
moves through the common enve l ope cr e ated by the ejected 
material (e.g. iMacDonaldl Il980l . Il986l ; IShara et alii 19861 ; 
iLivio et al.l fl99nir This causes the separation of the system 
to decrease and so the period decreases too. They found that 
the separation change given in equation becomes 



Aa Ami 
IT M 
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(14) 



l|Shara et al.lll986T ). where x = y/2 for a very slow nova and 
x = 1 for a fast nova. We note that ILivio. Govarie fc Ritterl 
(1991) suggest that assumptions made in deriving this mean 
that is is valid only for novae with longer orbital periods. 
Again, we can find the period change for this model with 
equation l|8|). In Fig. [1] we plot ^^i^m as a function of the 
mass ratio for this model for a slow nova (short-dashed line). 
The separation change (and hence period change) here is, as 
expected, less than that for the previous case without the 
frictional angular momentum losses. 



2.3 Magnetic Braking and Gravitational 
Radiation 

There are also two mechanisms that remove angular momen- 
tum from the system between outbursts that we consider 
briefly here. For systems with relatively long orbital periods 
(P > 3hr), magnetic braking provides the largest continual 
loss of angular momentum from the system. The rate of loss 
of angular momentum is given roughly by 



Jmb 
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(nyr) 3 M R|yr- 2 



(15) 



ijRappaport. Verbunt fc JossH l983). For reasonable parame- 
ters this gives a time scale of around a few billion years (e.g. 
iMartin fc Touill2005l ). 

For the closest systems (P < 3 hr) , gravitational radia- 
tion is the dominate cause of angular momentum loss from 
the binary system between outbursts. The rate of loss of an- 
gular momentum from two point masses in a circular orbit 
is 



JgR = 



32G i M1M2M 



(16) 



5c 5 a 4 

l|Landau fc Lifshitd[l95lh . This gives a timescale of order 
several million years except for the very closest of systems. 
Neither gravitational radiation nor magnetic braking will 
have any effect on a system on an observable timescale. 



3 COROTATING EJECTED MATERIAL 

In this Section we consider an entirely new mechanism that 
can change the orbital period during a nova outburst. Sup- 
pose ejected material that moves within a critical radius 
(that depends on the magnetic field strength) of the sec- 
ondary star couples to its magnetic field and so is forced to 
corotate with the binary orbit. The transfer of angular mo- 
mentum between the ejected material and the binary orbit 
causes a change to the orbital period of the system. 
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Figure 2. The change ^^^m as a function of mass ratio, q. 
For the solid line the ejected mass carries away only its specific 
angular momentum. All other lines take into account the angu- 
lar momentum change of the ejected material due to a magnetic 
field on the secondary (corresponding separation change in equa- 
tion [29} ■ In order of decreasing height the lines have R\ = 0.25 a 
(dotted line), R\ = 0.5 a (short-dashed line), R& = 0.75 a (long 
dashed line) and Ra = 0.9 a (dot-dashed line). 



The magnitude of the dipole magnetic field of the sec- 
ondary is 



B = 



(17) 



where [i is the dipole moment and R is the distance to the 
secondary star. The magnetic field energy density is 

B 2 

Emus = - — • (18) 
47T 



The kinetic energy density of the ejected matter is 

Emu = \pu\ (19) 

where p is the density of the ejected material. The ejection 
veloc ity, u, is in the range 300 - 3000 km s" 1 ((Sh ara et alj 
1986). From the continuity equation, close the the secondary 
star, we approximately have 



pu - 



M 
4ma? ' 



(20) 



where M is the mass ejection rate from the primary white 
dwarf star. We approximate the average mass-loss rate with 

Ami 



M 



(21) 



where r is the timescale over which the mass is lost. 

The magnetic energy density is equal to the kinetic en- 
ergy density at the Alfven radius which we find to be 



Ra = 



2/iV 
M11 



(22) 
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We consider this further in Section [4] where we apply our 
model to the recurrent nova U Sco. We parametrise the 
dipole moment with the magnetic field strength at the stellar 
surface, Bstar, so that 

(U = Bstar-Ra, (23) 

where R2 is the stellar radius of the companion which we 
find with equation (|13[) . 

The specific angular momentum of the ejected mass 
that is forced to corotate with the secondary is {a% + KR\) Q, 
assuming that the secondary is tidally locked. The distance 
from the centre of mass to the secondary star is 

Mi 



a 2 



M 



(24) 



The constant K depends on the distribution of the material 
within the Alfven radius. If it were distributed uniformly 
within a spherical shell, then K = 2/3. Because the density 
within the shell varies and the shell itself will not be perfectly 
spherical, we take K = 1. Thus, the angular momentum loss 
from the binary to the ejected material is given by 



AJ C , 



-/Ami(o 2 +Ra - o x )n. 



(25) 



We estimate the fraction of the ejected mass that gains an- 
gular momentum in this way to be 



f 



4a 2 ' 



Now with equations ((2)1 and (|25[) we find 



AJ C , 



J 



/Arm (1 + g) 2 
M q 



Ha 
a 



l + q 



(26) 



(27) 



We substitute this into equation ([6|) with 



(28) 

AMi = —Ami and AAI2 = to find the change in separa- 
tion 



Aa Arm J (1 + g) 2 

a M 



Ra\ 1-q 
a J l + q 



(29) 



We note that this reduces to equation ([7]) when there is 
no magnetic field. In Fig. [2] we plot ££{/' M as a function 
of the mass ratio for different values of the Alfven radius. 
The larger magnetic fields are even capable of producing an 
overall decrease to the orbital period during the nova for 
small mass ratios. 

In Fig. [l] we also plot SmWSf ^ or an Alfven radius of 
0.75 a (long dashed line) for comparison with the other mod- 
els we have considered. For this strong magnetic field, the 
period change is smaller than with mass transfer to the 
secondary or frictional angular momentum losses. If a sec- 
ondary star has a large magnetic field then it will signif- 
icantly alter the orbital period change during a nova out- 
burst. 

There is a critical mass ratio where the additional term 
in equation (|29[) causes no change to the separation (or or- 
bital period) 

(30) 



Qcrit 



i+(4r) s 



1 



If q < gcrit then the ejected material that is forced to coro- 
tate with the secondary gains angular momentum, angular 
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Figure 3. The solid line shows the critical mass ratio, g cr j t , as a 
function of the Alfven radius. If q < q cr ^t the period decreases due 
to the magnetic field, whereas if q > q cr jt the period increases. 
Below the dashed line, the period change becomes negative. 



momentum is lost from the orbit and so the separation de- 
creases. Similarly the orbital period change is smaller than 
previously predicted. However, if q > g cr i t , then the material 
forced to corotate loses angular momentum and the orbital 
period increases. We plot this in the solid line in Fig. [3] We 
note that most obse rved systems have a mass ratio q < 1 
jRitter fc Koibll2003h because there is a critical mass ratio, 
that depends on M2, above which the mass transfer becomes 
unstable (e.g. ISmith fc Vande Putti 2006). Hence in most 
systems the orbital period will decrease when the effects of 
a magnetic field are considered. 

We also plot the dashed line for the mass ratio below 
which the period change is negative. We see that for the 
larger magnetic fields it is possible that the orbital period 
may actually decrease for all mass ratios. This effect could 
be significant even for larger mass ratios. Fricti onal angular 
mom entum losses only dominate for q < 0.01 (|Shara et al.l 
1986) so this new mechanism potentially has a greater effect 
on the orbital period change. 



4 RECURRENT NOVAE 

There are now ten recurrent novae kn own in our galaxy (the 
tenth one was discovered last year; IPagnotta et alj 120091 ) 
and one system in the LMC. In this group, U Sco has the 
fastest decline rate of the light curve in past outbursts, 
and th e shortest recu rrence period (11 yr since the last out- 
burst, ISchaefer|[200ll) . It has outbursts recorded in 1863, 
1906 1917, 1936, 19 45, 1969, 1979, 1987, 199 9 (|Schmeerl 
Il999l ; ISchaefedlioiOal ) and 2010 (|Schaeferll2010bl ) and others 
have likely bee n missed because of its proximity to the Sun 
l|Schaeferll2004T ). 

U Sco has a white dwarf with a mass Mi = 1.55 ± 
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Figure 4. For a system with the parameters of U Sco. Left: The Alfven radius, R\/a, as a function of the magnetic field strength at the 
surface of the secondary, B s tar- Right: The mass fraction, /, that is ejected within the Alfven radius of the secondary. The upper three 
lines are for the upper limit on the nova timescale, r, and the lower three lines for the lower limit. In both plots, the solid lines are for 
ejection velocity u = 300 kms" 1 , the dotted lines for u = 1000 kms -1 and the dashed lines for u = 3000 kms — 1 . 



0.24 Mq (|Thoroughgood et all l200lh . We take the mass 
to be close to the upper limit for that of a white dwarf 
that is accretin g matter before a s upernova occurs, so 
Mi = 1.37 Mq (IHachisu et all l2000al ). This mass is con- 
sistent with the fact that U Sco has freq uent outbursts . 
The orbital period i s P = 1.2305521 d ijSchaeferl 1 19901 
ISchaefer fc Ringwaldl [l995h . The radius of a non-rotating 
white dwarf is given approximately by 



Ri = 7.99 x 10 8 



Mi 



Mi 
~Mch 



(31) 



where M c h = 1.44 Mq is the Chandrasekhar mass 
(|Nauenberg||l972l ). so the radius of the white dwarf in U 
Sco is Ri = 0.003 R . 

With equation {1} we find the mass accumulated before 
the outburst to be Ami = 2.36 x 10~ 6 Mq, consistent with 
estimates for the 1 999 outburst fkachisu et al.ll2000bf) an d 
the 2010 outburst (jBaneriee et alj|20ld : iDiaz et alj|2010h . 
We assume that all of this mass is ejected in the outburst. 
With this we find the ratio Ami/M = 1.05 x 10" 6 . The 
companion to the white dwarf in the system is a subgiant 
(|Schaeferi ri990) with a mass of M2 = 0-88 M and a radius 
of R 2 = 2.1 R (|Thoroughgood et al.ll200lT ). We can take an 
upper limit of r w 3 months (the timescal e on which the op- 
tical light curve drops back to quiescence. iMatsumoto et al.l 
I2003T I and a lower limit of r ~ a/u (the binary crossing 
time). We consider a range of ejection velocities of the ma- 
terial between u — 300 kms -1 and 3000 kms -1 . With this 
range of velocities, the lower limit to the ejection time is in 
the range 0.5 — 4hr. 

In the left hand side of Fig. U we plot the Alfven ra- 
dius (given in equation [22]) for this system for different ejec- 
tion velocities. For strong surface magnetic fields of a few 



thousand Gauss, it is a significant fraction of the binary 
separation. On the right hand side we show the fraction of 
mass that is ejected within the Alfven radius (given in equa- 
tion [26]) of the secondary star for different mass ejection ve- 
locities. Even for the strongest surface magnetic fields, the 
fraction of the ejected mass is always a small fraction of the 
total mass ejected. 

The mass ratio in U Sco is around q = 0.64. With Fig. [5] 
we see that with a large surface magnetic field, and hence 
Alfven radius of the secondary star, the period change dur- 
ing an outburst in U Sco is significantly reduced. With a 
strong enough magnetic field the period change could even 
be negative (see Fig. [3]). Measurements of the orbital period 
after the 2010 outburst are strongly encouraged in order to 
test this mechanism. 



5 DISCUSSION AND CONCLUSIONS 

We have considered the effect of a strong magnetic field on 
the surface of the secondary star on the orbital period change 
during a nova outburst. For most systems, the ejected mate- 
rial gains angular momentum as it couples to the magnetic 
field. The binary system loses angular momentum and the 
period change that results purely from mass loss decreases. 

Our results show that, contrary to expectations, the or- 
bital period in nova systems could decrease during outbursts, 
even in systems in which it would previously have been ex- 
pecte d to increase (e.g. forP > 8hr,[L ivio. Govarie fc Ritterl 
Il99lf ). Very accurate determinations of pre- and post- 
outburst orbital periods in recurrent nova systems are there- 
fore strongly encouraged. These could provide meaningful 
constraints of such quantities as the mass of the accreted 
envelope and the secondary's magnetic field. 
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Magnetic fields of the order of a few kilogauss 
on the secondary star have been suggested previousl y 
(e.g. iMever-Hofmeister. Vogt fe Meveri[l996l ; I Warner! [l996h . 
While such strong fields are typical of magnetic Ap stars, 
they may be less common in the secondaries of cata- 
clysmic variables. However, high magnetic fields have been 
discussed for cataclysm ic variables of shorter periods by 
iMeinties fe Jurual (2006). The subgiant companion in U Sco 
is expected to be synchronously rotating with the orbit, 
with a period of 30 hr, which is very fast for a subgiant. 
However, studies of late-type stars show that high fields can 
be ex pected for fast rotators (e.g. iNoves. Weiss fc Vaughanl 
1 1984 ) . 
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